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.Abstract 

This Quarterly Report deals with the deployment maneuver of a single-axis, 
vertical constellation with three masses. A new, easy to handle, computer code 
that simulates the two-dimensional dynamics of the constellation has been Imple- 
mented. This computer code is used for designing control laws for the deploy- 
ment maneuver that minimizes the acceleration level of the low-g platform during 


the maneuver . 
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1.0 I NTRODUCT I ON 

This is the first Quarterly Report submitted by SAO under contract NAS8- 
36606, "Analytical Investigation of the Dynamica of Tethered Constellations in 
Earth Orbit (Phase II), H Dr. Enrico Lorenzini, PI, and covers the period from 22 
February 1985 through 21 June 1985. 
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2.0 TECHNICAL ACTIVITY DURING REPORTING PERIOD AND PROGRAM STATUS 

2.1 Introductory Remarks 

The results of the study performed under NASA contract NAS8-35497 on the 
dynamics of tethered constellations in earth orbit were very encouraging for a 
single-axis vertical constellation. In the three-mass configuration the middle 
mass, if appropriately placed at the orbital center of the system, can be con- 
veniently utilized as a low-g experiment platform. A preliminary estimation of 
the acceleration level achievable was around 10~ a g when the system began its 
station-keeping phase from a rest condition and the J a term of the gravity field 
was taken into account. The above referenced report dealt with the station- 
keeping of the constellation only. No indications were provided on how to 
deploy the system to the final, steady state configuration. This Quarterly 
Report addresses problems related to deployment of a single-axis constellation 
with three masses from the initial, compact configuration to the final alignment 
of the constellation with the local vertical at fully deployed tether length. 


2.2 Two-Dimensional Equations Of Motion For A Single-Axis Vertical Constella- 
tion With Three Masses 

The SKYHOOK program was used to simulate the station-keeping phase of the 
constellation in Phase I studies. Station-keeping is easier to model than de- 
ployment so that a limited number of simulation runs were required. Deployment 
of a three-mass constellation is a new topic that requires a longer and more 
difficult effort. A complex computer code, like SKYHOOK, is therefore prohibi- 
tively cumbersome for a parametric study like this. For this reason a more 
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flexible and easy to handle mathematical model has been derived to study the 
deployment dynamics of a single-axis constellation with three masses. 

The model is two-dimensional . This provides a good approximation of the 
deployment maneuver that is dominated by forces acting in the plane of the 
orbit. The model assumes three point masses generlcally located in space (on 
the orbital plane) connected by two massless tethers of known elastic character- 
istics . 

The center of mass of the system is supposed to ha on a circular or nearly 
circular orbit but in any case the orbital motion is decoupled by the system 
oscillation. Presently no aerodynamic forces are in the model but they could be 
added easily. Aerodynamic forces, however, are not so important during deploy- 
ment at 500 km altitude. Excessive complication of the model would be impracti- 
cal since this model is developed in order to design a suitable deployment 
strategy for verification by SKYHOOK. 

The equations of motion were derived by using the Lagrangian formulation. 
The Lagrangian coordinates were chosen, as shown in Figure (2.2.1) , to give an 
immediate representation of the system oscillations. They also have the advan- 
tage of representing oscillatory modes that are, in some cases, very weakly 
coupled and therefore suitable for consistent simplifications in order to 
achieve analytical results. 
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2.2.1 General Case - 

The three masses are generically located on the plane of the orbit. With 
reference to Figure (2.2.1), the position of the masses with respect to a refer- 
ence frame centered at the center of mass and co-rotating with the system is 
given by: 


xi = (Rjfj - Ri*i + 1 1 ) sin 6 - R a « cos 8 


x a = (Rjf ] - Ri*i) sin 8 + (R x + Rj) < cos 8 


Xj = (Rjlj - Rif i — f j) sin 8 + Rj€ cos 8 


Z\ — (R jij — Rif i + fi) cos 8 + Rjt sin 8 


z a — (Rjf j — Rif i) cos 6 — (Ri + Rj) < sin 8 


z | = (Rjfj - Rifi - fj) cob} + R a « sin 8 (2.2.1) 


where : 


Ri — mi/m eo t 1 — 1 1 2 , 3 

Equations (2.2.1) are valid for t « f i( which is applicable in most cases. 
The potential energy of the system is given by: 
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t 

x 1 


Lagrangian coordinates ; 

0 = in-plane angle 

e = lateral defle-tion 

= tether length of tether a l a (orbit semi-major axis) 

= tether length of tether c 2 

«f to the center of the Earth 



Figure 2.2.1 Geometry and reference frame for a single-axis 
tethered constellation with three masses. 
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V 



n>i/ri 


( 2 . 2 . 2 ) 


which to the second order of approximation is as follows: 

V - - * £ a m t 

■ V 


1 2 » 

Zj _ z t + Xj - 3zj 

a 2a* 


(2.2.3) 


where (i is the gravitational constant of the earth and a is the semi-major axis 
of the orbit of the system center of mass. After considerable rearranging, the 
potential energy can be expressed as: 


= - | 0 J {m tot [2a J - 2(Rjf, - M i) 2 (3 cos J 0 -1)] + 


+ R a (m x + mj) [<* (3 sin J 0 - 1)] + 


+ mi[fj (3cos J 0-l) + 6fiRjt cos 6 sin 0 j + 


+ nj[fj (3 cos J 0 - 1) - 6t jRjt coc 0 sin 0 ]| 


(2.2.4) 


where (1 is the mean orbital rate. 


The kinetic energy of the system is given by: 

t = * m *M J 

* i 


(2.2.5) 


where : 


vi| J = [*! + 0 (a - z t )] 2 + [z t + Oxi] 2 


( 2 . 2 . 6 ) 
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After taking the derivative of equations (2.2.1) and aubstituting 
equations (2.2.6) and (2.2.5) we finally get: 


T * | {m t „[n J a J - (R,*, - R^j) 1 - (R,f, - R^)* (6 - 0)*] + 


+ Ra (®i + bji)[<* + O a « J + — 0) a ] + 


+ «i[^i + f! (4-n) a - 2(m 3 (/ lt - fio] + 

♦ mj[/j + t\ (i - ny + 2nRj(/,t - /,<)]} 


The Lagrangian function is easily computed as: 


L = T - V 


The equations of mution are given by: 



where the Lagrangian coordinates are: 


^ J 91a = « ; qj = i\\ = £ j 


The Qi's are the generalized forces given by: 



3£i 

3qj 


them into 


(2.2.7) 


( 2 . 2 . 6 ) 


(2.2.9) 


( 2 . 2 . 10 ) 


( 2 . 2 . 11 ) 


In equation (2.2.11) are the forces acting on the m 


and the radius 


tlu 
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vactora, respectively glvan by: 

?i ■ x t I + *1 k (2.2.12) 

Fi *= -Tx(coa^ + (/<ialn^)k - Tj (sin # - t/tico»0)l 

2 = [(Ti - Tj) coa 6 + (Ti/^i + T,/f,) c ain $ ]k + [(Ti-T,)sin 6 - (T x /t x + 7,/t t ) t cos 6 ]l 

* Tj(cosf - c/l|Sin0)k + Tj (sin 6 + «//, cos 0)1 (2.2.13) 

By substituting (2.2.12) and (2.2.13) into (2.2.11) we finally gat th* general- 
ized forcas as follows: 

Q* = o 

Q< = -<(T Jt x ♦ Tj/f j) 

= -Ti 

Q«i * -T, (2.2.14) 

whera Tj and Tj ara th a tensions in the respective tethers. 

By substituting the above into th* Lagrangian equations (2.2.9) and after 
performing the necessary derivations we finally get th* equations of motion as 
follows : 


<"i + (Rj*i - Ri*i)] + 2fi(4-fl)[^i + (Rj^i - Ri^i)] + 

+ 30 1 *! cos t sin 9[l x + 2 (R,f , - R x £j)] - cos (2#)} + 
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+ 21 ,( 4 -fl)[*,+ (R»*i -R^i)] + 
+ 30 1 ^i cos 6 aln 6[l , + 2 (K,t , - Ri^i)] + 3(J*f jR 2 c cos (20)} + 
■f Rj (m : + mj) [ c *§ + 2ti (6 - 0) - 3fl , t 1 cos 9 sin 6 ] = 0 


Rj(n>i + mj) [f - « (4 - tl) J - 30 J c sin J *] + 

♦ RjJHft[3n J co8«aintf (Rj/j-R^i) - 20 (Rj 4 j-R^) ] - *T,/i ,) (2.2.15) 


n>i{*i + (Rjij - Ri/i) - n J (3co» J « - l)[fj + 2(R,f, - R^,)] - 
- (4 - + (Rjfj - Ri^i) ] — 3ll , R J < cos 6 sin 6 - 2nRji| = - T, 


m,{/, - (Rj^i - Rx* x ) - 0 J (3cos^ - l)[f, - 2(R,f, - R» l,)] - 

- (4 - n) J [*, + (R,*, - Ri^o] + 30^2* cos e sin 6 +2nRj<} = - T, 

These equations of motion will be used to simulate the dynamics of the three- 
mass system during deployment as will bo shown later in this report. After 
further simplifications, relevant to particular cases, these equations can also 
provide information on the dynamics of the system without being numerically 


integrator! . 
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2.2.2 Case Of The Middle Mass At The System Center Of Mass - 

In equations (2.2.15) the term Rjfj - R x fi is the distance (along the line 
connecting mi and m 3 ) of the middle mass mj from the system center of mass. If 
the middle mass is maintained constantly at the center of mass by means of an 
appropriate control law then the above mentioned term and its derivatives can be 
eliminated in equation (2.2.15). The distances and lj can also be expressed 
as a function of the overall tether length t. From equation (2.2.13), by re- 
quiring that the component of Fj along the line connecting m x and mj is zero, we 
get T x = Tj = T so that equations (2.2.15) reduce to: 

m«, [f 1 6 + 2 It (0-O) J -I- 3Cl 2 t* cos 0 sin *] + 

+ R 2 (mi + raj; [« + 2€«(^-0) - 3fJ J « 3 cos 6 sin 0 ] = 0 

Rj[<-«(^-n) a -3O J esin J 0] = -eT/ffm*,) (2.2.16) 

m „[£ - 0 H (3 cos 3 0 - 1) - (6 - 0) J f ] - 
- Rj (mi - m 3 ) [3fl J e cos 8 sin 6 + 20 1 ] = — T/m.q 

Equations (2.2.16) provide some useful information on the dynamics of the sys- 
tem. The coupling among the rotational motion (equation in 6) and the lateral 
deflection (equation in t) is strongly reduced by placing the middle mass m a at 
the system center of mass. In particular if the initial lateral displacement is 
zero it stays zero It spite of the system libration. This suggests, for exam- 
ple, that the best way to deploy the constellation without overly perturbing the 
low-g platform Is to design a control law that keeps the middle mass at the 
system center of mass throughout the entire deployment. At the end of the 
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maneuver the mass mj can be moved to the orbital center (a few meters below the 
c.m.) in order to minimize the vertical acceleration. 

If we assume that the constellation is fully deployed, that 6 and 6 are 
equal to zero so that the tension is approximately constant and is given by T = 
3m aq n J £ , then the e equation in (2.2.16) reduces to: 

Rjl + n J e(3-R 2 ) = 0 (2.2.17) 

From this extremely simple equation we can derive an analytic formula to compute 
the free oscillation frequency of the lateral vibration of mass m 2 as follows: 

f« = -f- (3/Rj — l) 1/J n (2.2.18) 

As an example, for the constellation studied in Phase I, with: ni! = 90600. kg, 
m 2 = 4530 kg, = 9060. kg, O = 1.1068 x 10° rad/sec (altitude = 500 km) we 
get : 


f t = 1.45 hz ; T, = 688.4 sec (2.2.19) 

Notice also that in the first approximation the lateral oscillation frequency 
does not depend on the overall system length. This is because the tension is 
directly proportional to £ while the angle between the tethers and the line 
through the end masses mi and m 3 is inversely proportional to it. The tether 
length, therefore, cancels out when computing the horizontal component of the 


force acting upon the middle mass. 
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2.2 Numerical Integration Of The Two-Dimensional Equations Of Motion 

A new software code called 3MASS has been implemented to simulate the dy- 
namics of a 3-mass constellation by using the equations of motion previously 
described. The computer program uses a fourth order Runge-Kutta integration 
subroutine with variable step size. This integration has been successfully 
tested in the past with a similar mathematical model. Additional subroutines 
called DEPL0Y3 and STAKE3 model the tether control laws during the deployment 
and the station-keeping maneuver respectively. All the geometrical and mass 
properties and orbital characteristics of the constellation are read as input 
parameters by the main program and passed to the integrator. The integration 
subroutine, called DRKGS , calls the subroutine FCT which computes the deriva- 
tives of the Lagrangian coordinates. DEPL0Y3 or STAKE3 are called by FCT and 
after a successful integration step the output subroutine OUTP is called. OUTP 
formats the output variables for the printout and it prepares the output file 
for the plotter. Subroutine O'JTP also stops the program in the case of Drogram 
malfunctions or the achievement of the desired final conditions. 

Typical program malfunctions are: the in-plane angle is too large so that 
the constellation attitude is close to the horizontal plane or the integration 
step size becomes too snail. This computer program is very efficient, at least 
when the tethers are considered unstretchable . Normal maximum integration step 
sizes are around 50 sec. These characteristics make it suitable for deriving 
appropriate tether control laws for the deployment phase. 
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2.4 Deployment Maneuver 
2.4.1 Control Laws - 

Deriving a deployment control law depends on the control parameters avail- 
able to the deployer and on the measurement accuracy of these parameters. 
Knowledge of the in-plane and out-of-plane angular rate of the tether libration 
is required for a very effective damping of librations. Simpler control laws 
rely on the tether length and length rate or on length, length rate and tension. 
They are respectively called rate control laws and tension control laws even if 
the latter denomination is inaccurate. Rate control laws are easy to implement 
and are very reliable. However, they do not damp librations: the decrease of 
the libration amplitudes in this case is due to the conservation of angular 
momentum while the tether length is being varied. The tether tension is, on the 
contrary, related to libration angles and angular rates (well coupled to the in- 
plane oscillation and very weakly to the out-of-plane) so that a tension control 
law, if appropriately tuned to the tether longitudinal oscillation, provides a 
good damping of the in-plane oscillation. A major drawback is that the tension 
for short tether lengths is too low to be a reliable parameter to control. For 
this reason rate control laws (with or without angular feedback) have to be used 
for short tether lengths. In this study we used a rate control law for deploy- 
ment mainly for two reasons: a) if low level acceleration can be achieved in 
the middle mass, notwithstanding a residual oscillation of the system, then any 
other control law will be more effective than this control law; b) a rate 
control law constitutes the base for any other more refined control law. 

The system starts with the three masses very close together, though not 
coincident, orbiting in relative equilibrium. The lower mass, mj , is the space 
station, the central mass, mj, is the low gravity facility, and the higher mass, 
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mj, is a ballast mass. Thereafter the natural length of each cf the two tethers 
is derived from a control lav: 


i 


at ti < t < t T 

-0(t - h) l T < t < t at 


(2.4.1) 


where the various parameters may be different for the two tethers. Here l is 
the deployed natural length, t z an initial length, f T a transition length at 
which the reel begins decelerating, l ax the desired final length, and l t = t n ♦ 
o is the length which would result at infinite time, a being a small length 
chosen so that the deployment takes finite time but the cutoff is not too 
abrupt. The constant a = 0.75 O sin (2 B c ) is chosen by analogy to a single mass 
deployment where it would result in an attitude angle during deployment equal to 
6 C ; and 0 = af T /(f ( -f T ) is such that the deployment velocity is continuous 
through the transition. This control law may be solved explicitly to give the 
deployed natural length as a function of time. Assuming the initial time is t = 
0 , we get 


t (t) 



ti e°* 

l t - (t t - *t) exp [- 0 (t - t T )] 


0 < t < t T (2.4.2) 
t T < t < tsx 


where t T = ( 1 /a) ln(fj/fi) is the transition time and tsu = t T + (l/£) ln[(f t - 

*x)/<x)] is the turn off time. 

The same time constant (tea = 1/a) is used for the two tethers in order to keep 
the middle mass at the same location throughout the maneuver. The same consid- 
eration applies to the time constant for the deceleration phase (t e 0 = 


1/0) ■ 
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All the characteristic lengths, such as if, It, a are in the ratio Isici/lsm in 
order to provide simultaneous transitions in the two tethers. 


2.4.2 Simulation Results - 

Numerous simulation runs have been performed after implementing the previ- 
ously described control law in the new simulation program. Different values of 
parameters, such as a, f T , a and l s k have been tested to obtain a set of parame- 
ters that results in a well behaved deployment with a low level of acceleration 
of the middle platform. 

The dimensions and the masses used in these simulation runs are equal to 
the values used in the final report of Phase I. Namely they are: mi = 90600. 
kg, m 2 = 4530 kg, m 3 = 9060. kg, i SKl = 909.09 m, t s K 3 — 9090.9 m, orbital 
altitude = 500 km. Notice that this mathematical model does not separate the 
orbital center from the center of mass since the gravity field is approximated 
to the second order. This does not detract from the model since during deploy- 
ment the major perturbations on m 2 are due to the coupling of the lateral de- 
flection with the libration and since this coupling is not substantially af- 
fected by the third order term of the gravity field. The third order terms of 
the gravity field affect the vertical acceleration components so that a mass 
located at the system center of mass will actually experience a vertical accel- 
eration different from zero. The conclusion is that the results of these simu- 
lation runs are accurate for the horizontal acceleration component (the most 
important during deployment) while the vertical component does not show the bias 
dependent on the offset between the center of mass and the orbital center. 
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The best set of cortrol parameters came out to be: f T1 = 400. m, Ijj = 
4000. m, <Ti = 15 m, <r 2 = 150 m, 1/a = 1874.14 sec (that Implies 6 e = 20°). The 
deployment was completed in 5.1 hours while the simulation was stopped at 6 
hours in order to plot part of a station-keeping (with no tether control) . Two 
sets of figures follow that show results of a successful (and improvable) de- 
ployment maneuver. Figure ( 2 . 4.1 a,b) and Figure ( 2 . 4.2 a,b) show the tether 
length and tether speed for tether No. 1 and No. 2 respectively (tether No. 2 is 
called L3 to homogenize the notation in the equations of motion). Notice that 
these quantities, as mentioned before, are in the ratio l ski/^skj and notice 
also the small discontinuity of the tether speed at the end of the deployment 
due to the cut off of the exponential deceleration phase. The deployment is 
Initialized with l n = 2 m and £j 3 = 20 m. This condition can be actually 
realized with two deployment booms w... are extended at the beginning of the 
maneuver. These initial tether lengths affect the duration of the maneuver by 
expediting the initial phase of the deployment. Figure ( 2 . 4 . 3 ) shows the in- 
plane angle 9 of the constellation vs. time. Because of the choice of the 
Lagrangian coordinates (see Figure 2 . 2 . 1 ) Figure ( 2 . 4 . 3 ) must be read together 
with Figure ( 2 . 4 . 4 ), which plots the lateral deflection of the middle mass. The 
constellation's initial attitude is 20° consistent with the selected steady 
state inclination angle during the first phase of deployment. For this reason 
the in-plane angle is constant throughout the entire acceleration phase and 
decreases to the final, residual oscillation amplitude during the deceleration 
phase. The residual oscillation can be reduced to the steady-state value forced 
by the J 3 , as shown in the final report of Phase I studies, by means of a tether 
controlled station-keeping phase. It could also be reduced during the deploy- 
ment deceleration phase by using a tension control law. This second level 
improvement will be investigated in the next reporting period if desired by 


NASA. 
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Figure 2.4.1 Tether length vs. time 
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Figure 2.4.3 Constellation's in-plane angle vs. time. 
Initial value = 20°. 
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Figure 2.4.4 Middle mass’ lateral deflection vs. time. 
Initial value = 0. 
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Figure (2.4.4) shows that the lateral deflection t of the middle mass stays 
very small if the deployment is initialized with zero lateral deflection. This 
is because the coupling between lateral oscillation and libration is very weak 
when the middle platform is located at the system c.m. as mentioned before. Any 
other location of mass mj creates larger lateral deflection and consequently a 
higher acceleration level. Figure (2.4.5 a,b) show the respective trajectory of 
the lower and upper mass. Negative numbers mean trajectory outwards from the 
earth. Figure (2.4.6) shows the tension in tether 1 which is almost exactly 
equal to the one in tether 2. Figure (2.4.7) is the horizontal component of the 
acceleration of mass m 2 . The horizontal acceleration component is smaller than 
the accuracy of our computer in double precision (the middle mass initial lat- 
eral defleccion is equal to zero) and for this reason it must be considered only 
as an order of magnitude. Vertical acceleration components and acceleration 
modulus are not shown here because they are comparable to the horizontal compo- 
nent. Notice also that the vertical component has no bias due to the center of 
mass-orbital center offset as mentioned before. 

It is also interesting to see what the dynamic response and the accelera- 
tion level would be if the system starts the deployment with a small alignment 
error. Some simulation runs were, therefore, performed starting the deployment 
with different lateral deflections of the system. Results from a simulation run 
with an initial lateral deflection t = 0.05 m (2 inches) are shown in the 
following figures. The oscillation amplitude of the lateral deflection stays 
approximately the same throughout the entire simulation as shown in Figure 
(2.4.8). The damping of this type of oscillation deserves further study. Damp- 
ing terms were not included in the control law this time and it would be inter- 
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Figure 2.4.6 Tension in tether #1 vs. time. 
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Figure 2.4.7 Horizontal acceleration component of the middle 
mass vs. time for an initial lateral deflection equal 
to zero. 
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••ting to ••• how effective a tansion control law could ba in damping thia 
oacillation out. Figures (2.4.9), (2.4.10) and (2.4.11) ahow tha horizontal 
componant , vartical component and modulua of tha accalaration on the middla 
maaa . All of than are of the order of 10 * m/sec J (10 ’ g) but these value* aro 
strictly dependant on the initial lateral deflection and in any case they should 
ba fine-tuned (and poaaibly improved) with the use of a tension control law. 
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Figure 2.4.8 Lateral deflection of the middle mass vs. time. 
Initial value = 0.05 m. 



Figure 2.4.9 Horizontal acceleration component of the middle 
mass vs. time for an initial lateral deflection = 0.05 m 
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Figure 2.4.10 Vertical acceleration component of the midlie 

mass vs. time for an initial lateral deflection = 0.05 m. 



Figure 2.4.11 Acceleration modulus of the middle mass vs 
time for an initial lateral deflection = 0.05 m. 
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2.4.3 Concluding Remarks - 

The simple rate control law derived in this study proved to be effective in 
deploying a single-axis constellation with three masses. By using the same time 
constant for tether 1 and tether 2 and simultaneous transitions from one phase 
to another, the middle mass m 3 can be maintained at the center of mass of the 
system. This strategy reduces the acceleration level induced on the middle 
platform during deployment. When the deployment is over, the low-g platform can 
be moved to the orbital center of the system in order to reduce the vertical 
component of the acceleration. In the system under investigation in this Quar- 
terly Report the orbital center is only 1.2 m lower than the center of mass. It 
is also important to minimize the initial lateral deflection of mass mj with 
respect to the line through mj and mj in order to achieve a low acceleration 
level . 

Tension control laws should be investigated to reduce further the final 
libration and most of all to understand how effective these control laws are in 
abating the lateral oscillation of the system. 


3.0 PROBLEMS ENCOUNTERED DURING REPORTING PERIOD. 


None . 
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4.0 ACTIVITY PLANNED FOR THE NEXT REPORTING PERIOD 

During the next reporting period the analysis of the deployment phase of a 
single-axis constellation with three masses will be further developed. Tension 
control laws will be investigated. The effect of elastic tethers will be also 
considered. The tether control law that proves itself most effective for de- 
ployment will be verified by SKYHOOK or by a similar high fidelity model spe- 
cialized to simulate the dynamics of three-mass constellations. 

In addition to the above, the dynamics of the constellation when the middle 
mass travels along the tether in between the two end masses will be studied with 
the mathematical model derived during this reporting period. This model is 
perfectly suitable for investigating the situation where the middle platform 
travels along the tether. 


